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From the microscopic theory, we derive a number conserving quantum kinetic equation, valid for 
a dilute Bose gas at any temperature, in which the binary collisions between the quasi-particles are 
mediated by phonon-like excitations (called "condenson" ) . This different approach starts from the 
many-body Hamiltonian of a Boson gas and uses, in an appropriate way, the generalized random 
phase approximation. As a result, the collision term of the kinetic equation contains higher order 
contributions in the expansion in the interaction parameter. This different expansion shows up that 
a scattering involves the emission and the absorption of a phonon-like excitation. The major interest 
of this particular mechanism is that, in a regime where the condensate is stable, the collision process 
between condensed and non condensed particles is totally blocked due to a total annihilation of the 
mutual interaction potential induced by the condensate itself. As a consequence, the condensate is 
not constrained to relax and can be superfluid. Furthermore, a Boltzmann-like H-theorem for the 
entropy exists for this equation and allows to distinguish between dissipative and non dissipative 
phenomena (like vortices). We also illustrate the analogy between this approach and the kinetic 
theory for a plasma, in which the excitations correspond precisely to a plasmon. Finally, we show the 
equivalence of this theory with the non-number conserving Bogoliubov theory at zero temperature. 

PACS numbers: 03.75.Hh,03.75.Kk,05.30.-d 
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I. INTRODUCTION 
A. Superfluidity and H-theorem 

A lot of studies have been devoted to the theoretical 
understanding of statistical and dynamical properties of 
a weakly interacting Bose condensed gas. In particular, 
many works have been accomplished on the derivation of 
quantum kinetic equations (QKE) that govern the evolu- 
tion of the condensate fraction together with his thermal 
excitations Q, S H SI S Si, E El. In a sim- 
ple microscopic model of an homogeneous gas, one can 
describe the condensate by the atoms that populate the 
lowest ground state of energy and the thermal excitations 
by the atoms contained in the excited energy levels. The 
population of atoms in each level evolves according to 
the probability of scattering between the atoms. In the 
case of a uniform gas, the so-called Uehling-Uhlenbeck 
quantum kinetic equation (UUQKE) is a non-linear in- 
tegral equation which describes the detailed balance of 
the population transfer of atoms for each mode of the 
wave- vector through a binary collision term [T^ . This 
term depends linearly on the scattering differential cross 
section and nonlinearly on the mode population. It has 
also the remarkable properties to allow the QKE to obey 
conservation laws. These laws guarantee that the total 
number, the total momentum and the total kinetic energy 
of particles are preserved during the collision processes. 
More striking is the law stating that the production of 
entropy must be always positive, guaranteeing that the 
system obeys the second law of thermodynamics and, 
consequently, that it is always dissipative. This impor- 
tant requirement is known as the Boltzmann H-theorem 
established for a classical gas. 

The UUQKE has been derived in a weak coupling ap- 



proximation, valid for a diluted gas for which the kinetic 
energy is much higher than the potential energy. The 
collision term is indeed a second order expansion in the 
interaction potential between atoms leading to an expres- 
sion of the differential cross section in the Born approxi- 
mation. 

One could use such a QKE in a regime below the crit- 
ical point of condensation. In particular, for an inho- 
mogeneous gas in a trap potential, the kinetic equation 
describes the evolution of the Wigner function and must 
contain additional terms, taking into account the free 
propagation of the atoms, the influence of the external 
potential and the Hartree-Fock mean field. The conser- 
vation laws are still valid and allow to express the hydro- 
dynamic equations, including the equation for the local 
entropy production. 

Despite these consistencies, the resulting QKE suffers 
from the lack of understanding of an important phe- 
nomenon: superfluidity If this phenomenon really 
exists for a dilute gas then the second order the- 
ory must necessarily be revisited as it does not take 
into account the frictionless motion. Indeed according 
to the H-theorem, the homogeneous gas with a zero to- 
tal momentum evolves towards a statistical equilibrium 
state characterized by the Bose-Einstein particle num- 
ber distribution — l/[exp(/3(ek — A*)) — 1] where k 
is the wave- vector, P = l/ksT the inverse temperature, 
^ the chemical potential, and ek = k^/2m the particle 
kinetic energy respectively (fi, = 1). For /i ^ and 
k — > 0, the ground state has a macroscopic population 
no — l/[exp(— /3/i) — 1] having no relative velocity with 
the non condensed part of the gas. A non zero relative 
velocity corresponds to a non equilibrium situation and 
collisions between condensed and non condensed atoms 
will irrevocably damp this velocity towards zero. 

From this observation, we conclude that the second or- 



2 



der theory is no longer valid as far as superfluidity is con- 
cerned. Obviously, at zero temperature, only the poten- 
tial energy is the dominant contribution since the conden- 
sate is at rest and no thermal excitation subsists. More- 
over, other indications confirm that an improved descrip- 
tion of a weakly interacting Bose gas requires a higher or- 
der analysis in the interaction parameter. Among them, 
let us mention that the entropy of the condensate must 
be zero or close to zero while the H-theorem predicts an 
entropy 5*0 ~ log no corresponding to a system in the 
Grand Canonical ensemble with large statistical fluctu- 
ations of the particle number compared to his average 
value 5no/no ~ 1 0. This is due to the Bose enhance- 
ment factor which stimulates the collision rate in a huge 
manner as long as a condensed particle is involved in the 
process. On the other hand, the equilibrium statistical 
approach based on the partition function formalism in- 
dicates that the presence of a small interaction lowers 
considerably the fluctuations to an amount irrelevant in 
the thermodynamic limit Sno/no — » |l6| . 

Landau was the first to give an explanation of the su- 
perfluidity mechanism. He found a necessary but not 
sufficient condition for which this phenomenon happens 
0, 0| . He showed that a superfluid interacting with an 
external body (for example the wall) cannot release its 
energy, unless it evolves with a velocity higher than a crit- 
ical one. The argument is based on the impossibility to 
satisfy the momentum-energy conservation requirement 
because the excitation emitted by the superfluid has a 
phonon-like dispersion relation. 



B. Beyond the second order perturbation theory 

Attempts to analyze contribution coming from higher 
terms have been carried out with some success ( see 
[H 113 for a review). At zero temperature, Bogoliubov 
has calculated corrections to the ground state energy. 
The result is non analytic in the interaction parameter, 
due to infrared divergencies, and is obtained through 
a re-summation of an infinite number of contributions. 
Moreover, the theory predicts the existence of a phonon- 
like excitation, necessary in order to justify the Landau 
mechanism for superfluidity. 

Using Green function techniques, Beliaev improved the 
description of the phonon like excitation namely by calcu- 
lating its damping rate. Later on, Hugenholtz and Pines 
(HP) demonstrated that this excitation must necessarily 
be gapless. If at zero temperature, the theory seems well 
understood and widely accepted, for flnite temperature it 
seems rather controversial. Popov was one of the first to 
extent the Bogoliubov and Beliaev works using the Mat- 
subara formalism at finite temperature 20] . However, 
this approach is essentially valid for a weakly depleted 
gas and thus for temperatures much below the critical one 
Tc- To take into account a strong depletion, Girardeau 
derived the Hartree-Fock-Bogoliubov (HFB) mean field 
equations which unfortunately, have, the strong inconve- 



nient of having a gap in the quasi-particle energy spec- 
trum PH . 

Since then, many attempts, have been made, in order 
to suppress this gap and "rescue" the HFB equations, 
with the purpose to go beyond the Bogoliubov phonon- 
like dispersion relation. [23 . l23l |23 | . A very popular one 
is the so-called " Popov approximation" which consists in 
suppressin g th e product of anomalous term in the HFB 
equations |22| . Another consists in the renormalization 
of the HFB equations but at the price of making the ap- 
proximation of a weakly depleted Bose gas jl^. Conse- 
quently, these approaches are only strictly valid for tem- 
peratures much below Tc- A derivation of a well-defined 
theory, which is able to describe the weakly interacting 
Bose gas for any regime of temperature, is still an open 
problem and remains to be established (see the conclusion 
of 25]). By well-defined, we mean that the expansion in 
the small interaction parameter must be valid, whether 
the condensate is strongly depleted or not. 

Nevertheless, using these approaches, QKE have been 
derived taking into account these higher order effects in 
the collision term. While some authors use the dispersion 
relation resulting from the Bogoliubov theory 6, 7, 26], 
others use the one resulting form the renormalized HFB 
theory Superfluidity can be achieved in these models. 
Indeed, the resulting QKE can preserve detailed balance 
even for metastable states, in which a non zero relative 
velocity persists. In this situation, the Bogoliubov dis- 
persion relation becomes asymmetric - due to a Lagrange 
multiplier representing the relative velocity - and remains 
non negative, as long as this velocity does not exceed the 
sound velocity. Otherwise, the dispersion relation is neg- 
ative and the gas becomes unstable jl^. Although these 
models represent a considerable amount of work, they do 
not yet provide a full account of the conservation laws 
and the H-theorem that might be deduced from their 
QKE. 

However, let us mention that a QKE has been proposed 
valid in principle for temperatures close to Tc 01 Ell- The 
main problem of this approach is that the dispersion rela- 
tion has a gap and thus fails to explain the superfluidity. 

In general, in all these models, collisions between con- 
densed and non condensed atoms are possible and their 
rate is huge because of the Bose enhancement factor in 
the same way as in the UUQKE. If we follow the same 
reasoning as above and if an H-theorem really exists, we 
might expect a condensate entropy of the same order of 
magnitude So ~ log no. As a consequence, the parti- 
cle number fluctuations of the condensed mode are also 
huge, in contradiction with estimations made from the 
partition function formalism 



C. The random phase approximation approach 

The present paper is devoted to provide a possible al- 
ternative explanation to the problem of superfluidity in 
a diluted Bose gas. For this purpose a different QKE 
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is derived in the random phase approximation (RPA) in 
which higher order terms in the interaction parameter are 
retained in the binary coUision term (in some work, the 
RPA QKE refers to the UUQKE and thus has a differ- 
ent meaning from this paper [5|). Furthermore, in com- 
parison with other approaches, the QKE is valid for any 
regime of temperature below and above since, under no 
circumstances, the approximation of weak depletion has 
been used. 

The RPA is commonly used to describe the collective 
excitations in the quantum plasma of an electron liquid 
■ The idea behind this approximation is to neglect 
contribution containing average over pair of field opera- 
tors that are not oscillating in phase but rather randomly. 
In the terminology of optics, we neglect contribution that 
are not phase matching. These contributions come essen- 
tially from averages over pairs of different modes which 
oscillate with a relative random phase. The reason for us- 
ing the RPA in a diluted plasma is that, we expect that 
the excitations propagate over a sufficiently long time 
that non phase-matching terms are destroyed by inter- 
ference. By analogy with optics, we make somehow the 
far field approximation. The analysis of collective excita- 
tions in plasma using RPA reveals that the coulombian 
interaction potential between the electrons is shielded at 
long distance, due to the dynamic dielectric function. In 
fact, these results generalize the Debye theory predict- 
ing the screening of the potential between a charged ion 
due to the presence of other surrounding ions. As an ex- 
tension, the RPA includes also the dynamical aspects of 
these charged particles. 

A QKE for the plasma has been derived '28', "^gll. It 
predicts that the coulombian interaction potential for the 
cross section is shielded precisely by this dynamic dielec- 
tric factor. The classical version of this equation is known 
as the BGL kinetic equation (Balescu-Guernsey-Lcnard) 
|l2j| . For deriving it, Balescu uses a re-summation of 
the so-called ring diagrams, which is another approach 
equivalent to the RPA 28]. Later on, Wyld and Pines 
established a connection between the QKE and the plas- 
mon theory j2p|. In their approach, the shielded poten- 
tial results from a more subtile dynamical mechanism, 
in which the two electrons must emit and reabsorb an 
intermediate plasmon during their interaction. The dis- 
persion relation of the plasmon corresponds precisely to 
the collective mode of the plasma and its decay rate to 
the Landau damping. 

In this paper, an identical type of approach will be re- 
produced in the case of a dilute Bose gas. It is already 
known that the RPA allows to recover the Bogoliubov re- 
sults for t he g round state energy but with one important 
difference Namely, in the RPA, there is no need of a 
spontaneous breaking of the symmetry U{1) but rather, 
during the derivation, the total particle number is kept 
conserved. Note that a number-conserving formalism has 
been already used to derive a QKE H- But this formal- 
ism is based on a expansion [33 and thus differs 
from the one in the RPA \m 



These considerations allow to derive a number con- 
serving QKE in which the interaction potential is also 
modified by the presence of a dynamic dielectric func- 
tion resulting from collective excitations. However, in 
contrast to a plasma, this function presents a strange be- 
havior according to the kind of quasi-particles interaction 
it mediates. When the interaction involves a particle of 
the macroscopic condensate, it has the effect to totally 
annihilate the potential, forbidding any binary collision 
to occur. 

This surprising result allows to explain why a 
metastable condensate cannot decay into a state of lower 
energy. Simply, it is forbidden for a particle of the con- 
densate to scatter with the excited particle and to become 
thus excited. No exchange of particle can occur between 
the two fiuids. This unexpected phenomenon is a conse- 
quence of the collective behavior induced by the presence 
of a macroscopic condensate. In short, collision cannot 
happen because an induced mean field force, generated 
by the condensate, compensates exactly the interaction 
force felt by the condensed and non condensed particle. 
The net result is the absence of an effective interaction 
force preventing any scattering. In other words, the di- 
electric function used to attenuate the binary interaction 
potential becomes simply infinite. 

This phenomenon of "collision blockade" also influ- 
ences the mechanism of condensate formation. Since par- 
ticle cannot be exchanged with the macroscopic conden- 
sate, other mechanisms must be found. A closer analysis 
reveals that the collision blockade happens only when the 
Bose gas is stable. By stable we mean that the collective 
oscillations are always damped. Precisely, the Landau 
damping plays this role in both a plasma and a con- 
densate as long as we are close to equilibrium 0, |^ . 
However, for some non equilibrium situation, it has been 
shown for a plasma that the oscillations can grow expo- 
nentially leading to an instability [s^lssf . Such behavior 
happens also for a Bose gas which becomes thus unsta- 
ble. For example, we will show that this is the case when 
the relative velocity between the normal and superfluid is 
higher than the speed of sound. In such instable regime, 
the picture of a collision blockade is no longer valid and 
the QKE gets more complicated. Such more sophisti- 
cated QKE can describe the exchange of particle with the 
condensate [s^, 113 and provides an alternative explana- 
tion for the condensate formation [s^ . l35l l3^ . Another 
possibility is to exchange them on the edge region where 
the condensate is not macroscopically populated so that 
scattering can occur. But this requires also a really spe- 
cific analysis. Therefore, the present QKE derived here 
will not address the important issue of particle exchange 
with the condensate. 

Assuming that collision processes are local in space, 
the derivation can be extended for a weakly inhomoge- 
neous Bose gas. In this way, we recover the generalized 
Gross-Pitaevskii equation for the condensate, but with 
the difference that the absence of a binary collision term 
in this equation has now a clear justification. When the 
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normal cloud is in a local thermal equilibrium, no dissi- 
pation exists anymore and we can derive a set of coupled 
equations defining the superfluid regime for any finite 
temperature below Tc. 

The paper is divided as follows. We first begin by 
a heuristic approach of the "collision blockade" mecha- 
nism in section 2. In section 3, we derive the QKE for 
a homogeneous Bose gas. We review the simple RPA in 
which only the Hartree or direct term is considered and 
the generalized RPA (GRPA) in which both Hartree and 
exchange Fock terms are considered An equation 

of motion for the particle-hole pair operator is derived 
from which the collective and scattering excitation energy 
spectra are deduced. The solution of the equation of mo- 
tion allows to calculate the collision term and to predict 
the collision blockade phenomenon. We show that the 
QKE conserves the total particle number, the total mo- 
mentum and the total energy (kinetic and Hartree-Fock) 
and that a Boltzmann-like H-theorem exists. In section 
4, the reasoning of the previous section is extended to 
the case of a weakly inhomogeneous Bose gas. From the 
resulting QKE, we delimitate the superfluid regime i.e. 
the regime in which there is no local production of en- 
tropy. In section 5, by analogy with plasmon theory, we 
interpret the collision phenomenon as the result of the 
exchange of an intermediate boson (condenson). In sec- 
tion 6, we show how the number conserving RPA theory 
allows to recover the Bogoliubov results for the ground 
state energy. Finally, section 7 is devoted to conclusions 
and perspectives 



II. THE "COLLISION BLOCKADE" 
PHENOMENON 

A. Preliminary definitions 



to renormalize the theory by removing eventual high en- 
ergy divergencies. In what follows we will concentrate on 
a repulsive contact interaction a > 0. 

The so-called annihilation and creation operators of 
the particle-hole pair are of interest: 



(3) 



They represent an excitation of momentum q created 
from a particle which transfers its momentum from k 
to k + q. The kinetic energy transfer for this excitation 
is given by Wk,q = Ck-i-q — Ck- In particular, we define 
the density fluctuation operator pq = Sk/^k.q — pLqj 
the number operator rik = /3k,o and the total number 
N = X^k^-k- In terms of these operators, the Hamilto- 
nian becomes: 



H 



k^ 

k 



2m 



N) 



(4) 



B. Heuristic approach 

Assume a macroscopic condensate containing rik^ par- 
ticle evolving with momentum kg. In the absence of 
excitations, there are no fluctuations of the density op- 
erator i.e. (/Oq)'^^ = (Jq.o'^ks- Suppose that an exter- 
nal potential is turned on creating locally fluctuations 
of the density of the condensate. Expressing these per- 
turbations in Fourier space, we can characterize the ex- 
ternal potential (/)ea;t(q, and the density fluctuations 
(5n(q, w) = (Spq) = (pq) — (Pq)'^'' in terms of its wave- 
vector q and frequencies lu components. The Hamil- 
tonian created by this external potential is given by 
Hext{t) = X;qPlje"'"'0ea;t(q,t^)+c.c. The linear rcspousc 
to this potential is given by the formula: 



We start from the Hamiltonian: 



^ = 5l«k4ck+ 5Z ^4+qCL-qCkCk' (1) 



k,k',q 



cjj. and Ck are the creation and annihilation opera- 
tors obeying the commutation relations [ck, cj^.,] 



^'k.k' 



[Ck,Ck'J = [Ck, cjj.,] = 0. Ck = 5^ is the kinetic energy of 
the particle. Uq is the interaction potential expressed in 
Fourier transform. Since we are concerned only with low 
energy binary collisions in the channel Z — 0, it is com- 
mon to replace the potential Uq by a pseudo-potential or 
contact potential that we treat in the Born approxima- 
tion |l3: 



Ana , Ana v-^ 

TYl Tn 1/ ' ^ 



1 



mV ^ ek+q + Ck'-q - Ek - Ck' 



(2) 



where a is the scattering length. Usually we consider 
only the flrst linear term in the scattering length. The 
second term is ultra-violet divergent and is only present 



(5n(q,Lj) = x(q,^)'/'ea:t(q,t^) 



(5) 



where x(q, is the susceptibility function. This re- 
sponse function can be calculated formally by treating 
Hext{t) as a perturbation in the first order. Without the 
presence of a binary interaction potential between the 
particle, this function is given by [27|: 



xo(q,w) 



(6) 



It represents the transition amplitude for a condensed 
particle to increase its kinetic energy by an amount 
i^ks.q = Cks+q — Eks minus the amplitude for an excited 
particle to be transferred to the condensate releasing the 
energy — Wks-q,q = ^k^-q — Ck^- An infinitesimal quan- 
tity iO+ has been added to ensure the convergence and 
is due to the adiabatic switching process of external per- 
turbation. In the presence of the interaction, the conden- 
sate acquires a self-interaction potential energy given by 
Uqhi^^I {2V). Moreover, according to the RPA, the pres- 
ence of fluctuation density changes the potential energy 
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by the amount SHint = X]q PqS '"'^0(^1 ^) + c-c. where 
(5(/)(q, cj) = (C^q/^^)(^Pq) is the potential induced by the 
presence of the density fluctuations. Consequently, the 
global response of the system becomes: 

Sn{q,u}) = xo(q,w) (0ext(q,w) +(50(q,w)) 

= Xo(q, ^) (0ext(q, ^) + ({/q/F)5n(q, w))(7) 

Comparison between (jSJ and (O allows to deduce: 



x(q,t^) 



xo(q,w) 



l-(C^q/V^)xo(q,c^) 



(8) 



The total potential created inside the system (j)tot (q, ^) = 
^ea;t(q, + S(f){q,uj) is related to the external potential 
through the dynamical dielectric function: 



/C(q,c.) 



/-tot q,w) y 



(9) 



In the electromagnetism language, the gradient of 
4>exti(lj i-^) corresponds to the electrical displacement, 
while the gradient of (ptot (q, (jj) corresponds to the electric 
field. These are respectively the external and total field 
force acting on the particle. The dielectric function usu- 
ally has the effect of attenuating the external field force 
by means of an induced force so that the total field force is 
smoothed out. In particular, we notice that this function 
is infinitely resonant for frequencies uj = Lo^^^q, — ij-'ks-q.q 
which means that the induced potential exactly compen- 
sates the external one. For these frequencies, the external 
potential affects the density fluctuations but has no lo- 
cal influence anymore on the condensate particle itself. 
This important observation is at the origin of the colli- 
sion blockade phenomenon. 

Indeed, assume that the source of this finite external 
potential results from the transition of an excited parti- 
cle of momentum k' towards another excited state with a 
momentum k' — q. The released energy during this pro- 
cess is cj = — Wk',-q- According to the Fermi golden rule, 
a collision occurs between the condensed and non con- 
densed ingoing particles if the transfer energy are equal 
— Wk'.-q = f^ks.q or equivalently the total energy is con- 
served +ek' = ^ks+q + ^k'-q- But then, the dielectric 
function becomes infinitely resonant A^(q, w) —>■ oo caus- 
ing the total potential 0tot(q, i^) — > 0. Thus the inter- 
action potential felt by the condensate particle becomes 
non-existent, since the induced mean field cancels the ex- 
ternal potential generated by the excited particle. That 
precise phenomenon is responsible for the absence of an 
effective scattering process since, without interaction po- 
tential, no scattering amplitude can appear. The same 
effect happens when the scattering involves an outgoing 
particle in the condensate whereas, in such a case, the 
transfer energy — Wks-q,q causes the infinite resonance. 
Thus, these giant resonances have the effect to protect 
the condensate particle from scattering with the others. 
We must notice that this mechanism works only if the 
condensate population is macroscopic, otherwise, no in- 
duced potential is generated in the thermodynamic limit 



since, for rik^/V 0, /C(q, w) — > 1. As we shall see 
in the next sections, a more elaborated model confirms 
this prediction in the more general case of a non equilib- 
rium Bose gas, in which the condensate can be strongly 
depleted. 



III. KINETIC THEORY IN THE RPA 

Many methods have been developed to derive kinetic 
equations for a dilute Bose gas ilj, . In order to arrive 
rapidly to the final result , instead of using complicated 
many body techniques, we base the derivation on an op- 
eratorial method developed in [23| which appeared to be 
a simpler way to reach results without loss of generality. 
The approximations are the foUowings: 1) homogeneous 
Bose gas, 2) thermodynamic limit, 3) generalized RPA, 
4) instantaneous collisions (Markovian QKE), 5) no frag- 
mentation of the condensate (only one macroscopic pop- 
ulation mode). 



A. The random phase approximation 

In this subsection, we give a brief overview of the RPA 
developed in p^ . This approximation has been used 
quite extensively for the quantum electron liquid for de- 
scribing the screening effect. 

The dynamic of pk.q is given by the Heisenberg equa- 
tion motion. Using the relation: 

[Pk,q, Pk%q'] (^k+q^k'Pk^k'+q'-k — '5k-q',k'Pk',k+q-k(10) 

we find 



i^Pk,q = [Pk,q, H] = (ek+q - ek)Pk,q 
+ ^[/'k.q+q' - Pk-q',q+q',Pq,] + 



q' 



(11) 



where the brackets refer to to an anti-commutator. Two 
cases in the RPA are generally considered in systems close 
to equilibrium j^^. The simple RPA, in which only the 
Hartree or direct terms contribute, and the generalized 
RPA, in which the Fock or exchange terms are also re- 
tained. We shall concentrate on the second approxima- 
tion, since for a contact potential Hartree and Fock terms 
are identical. Only for the condensate-condensate inter- 
action, the Fock term does not appear. In the GRPA, 
we consider that operator pk.q with non zero momentum 
transfer q 7^ gives a negligible contribution in compar- 
ison to fiki as it involves different modes oscillating with 
a random relative phase. 

The procedure to get the non equilibrium RPA equa- 
tions is as follows. For a momentum transfer q = the 
Ea. (|ll|) is kept unchanged. For q ^ 0, however, we keep 
among all terms those combinations of creation and anni- 
hilation operators involving product of Pk',q and nk" for 
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all possible values of k' and k", and neglect those com- 
binations that cannot be written in this form. These re- 
moved contributions are quadratic in the operator pk'.q' 
with q' ^ q, 0. In this approximation, only contributions 
conserving the momentum transfer are relevant, the oth- 
ers coupling the various pk,q with different momentum 
transfers are neglected. We expect that the excitations 
of momentum q propagate, without interfering with the 
others, over an enough long time determined by the di- 
lution of the gas, that (this would correspond to the far 
field limit). The result is for q = 



E 

q' 



2V 



k,q' — Pk- 



(12) 



and for q 
. d 



(ek+q - ek)Pk,q 

V-^rC^k-k'+q , 



k' 



2V 



k-k' 



^ 2V 

k^k 



Fk 



E 

k'#k-q 

E 



f/k- 



2V 
2V 



2V 

»^k+q, Pk',q] + 
[«k,Pk',q] + 
[^k+q,Pk',q] + 



(13) 



Ea. (|13|l is an integral operatorial equation linear in pk.q- 
We can linearize this equation by averaging all possible 
bilinear contributions. Since the homogeneity of the gas 
imposes (pk.q) = if q 7^ 0, we are left with the average 
on the number operator (rik) — 7^k and we obtain an 
integral equation which possesses the same structure as 
Eq.(5.183) p. 318 of 0: 



d_ 

'dt 



i— - Cf^^^ 



Pk.q = ("k 



r^k+q) XI TTPk'.q 



k'T^k 



-rik 



E ■ 

k'7^k-q 



C/k-l 



V 



-Pk', 



k'5^k-|-q 



t/k-t 

V 



appear, we should be cautious not to count twice terms 
between the same modes. We have avoided this difficulty 
by excluding carefully in the sum over the modes those 
leading to a double counting and which appear only in 
the integral terms. Finally, another difference is that rik 
is still a function depending on the time, generalizing in 
this way the RPA approach for systems in non equilib- 
rium. 

For comparison, Zaremba et al. have derived a QKE 
in an approximation which would correspond to a differ- 
ent equation for Pk,q- As we shall see below, their anal- 
ysis is equivalent to removing the integral terms in l|14|) 
and to approximating the quasi-particle energy by 



^ZNG 



k^ 

2to 



1^ — "k' 



V 



E 



Uk-l 

V 



-Uk' 



(16) 



We notice immediately that they have avoided a double 
counting in the exchange term. As a consequence, in the 
particular case of a contact potential and for kg = 0, the 
energy difference that they obtained has a gap: 



^ZNG 
^0 



_ k^ Ana no k^o 4™ no , , 



2m 



m V 



m V 



(14) 



Indeed, in the term containing the bracket, we recognize 
the difference ujk,q — ^k+q ~ ^k^"^ between the quasi- 
particle energy calculated in Hartree-Fock approximation 
(HFA): 



Clearly, for k ^ and if the mode kg = is macroscopi- 
cally populated, this finite gap cannot be neglected in the 
scattering energy spectrum. In our approach, however, 
this gap does not appear anymore in £k^^- 

The two operatorial equations (|12l) and H14() are the 
equations leading to the QKE. 



B. Collective and scattering excitations 

It is instructive to analyze the frequency spectrum so- 
lution of the eigenvalue problem given by Eq. H14|) . In an 
electron liquid, the linear equation possesses two kinds 
of eigenvectors: (i) the scattering solutions involving the 
presence of only one particle and (ii) the collective solu- 
tions involving the presence of many particles. 

For reasons we will explain below, let us concentrate on 
the problem of calculating the impulsion response or di- 
electric propagator i^q(k, ki, t) to an initial particle hav- 
ing a momentum ki. We replace the operator pk,q by 
this function and we substitute the interaction potential 
by its first order expression in (|2J). After simplifications, 
the integral equation for this function is: 



HFA 



k^ 

2m 



E-^"k' 

k' 



EC^k-k' 
-^^rik' 

k' 



V 



(15) 



m 2m'^'' 

— [ V iYq(k,kl,t) 



k^k 



On the other hand, in the integral terms, both Hartree 
and Fock terms are present. However, Ea. H14() shows 
some differences: firstly, since we are dealing with bosons, 
the Fock terms have the opposite sign compared to 
fermions; secondly, since a macroscopic condensate might 



Y: Wq(k',kl,t)-!^ J2 Wq(k',kl,i)|l8) 



k^k-q 

with the initial condition: 



k'#k+q 



iYq(k, kl,i = 0)=5k,ki 



(19) 
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If nk varies slowly in time and thus can be consid- 
ered as constant during the evolution of the impulsion 
response, Ea. lfT^ can be solved exactly in the thermody- 
namic limit. Similar equations have been solved in the 
field of plasma physics ,4Qj. For this purpose, we define 
the Laplace transform as: 

/>oo 

Wq(k,ki,cj) = / dte*("+*+)*^^q(k,ki,i) (20) 
Jo 

A 0+ has been added in order to ensure convergence of 
the integral. If we suppose that kg is the wave-vector 
for the superfluid mode, then we can make the decom- 
position between a normal component and components 
involving the superfluid mode: 

Wq(k,ki,a;) = Wq(k, ki,tj) + 5k.ksWq(ks, ki, w) 

+'5k,ks-qWq(ks - q,ki,cj) (21) 

Also, we distinguish the normal mode population n'-^ = 
(1 — ^k,ks)"-k from the condensed mode. Plugging this de- 
composition into (|I4|) and neglecting some n'-^ by taking 
the thermodynamic limit, we obtain for the superfluid 
modes: 

[uj + iO+ - - ^]Uq{ks,ki,uj) = i(5ks,ki 

TO ZTO 



where 



(^2Z^q(k',kl,C^)+iYq(k„kl,0.) 



+Wq(ks-q,ki,c^)) (22) 



[ui + iO+ 



ks-q , q 



TO 2to 
7ra rik. 
TO V 



— ]ZYq(ks - q, ki , w) = i^ks-q,ki 



Ana n ' 

- — 17^(E 2Ziq(k',kl,C^)+Z^q(k„kl,C^) 



hi^q(kg-q,ki,c^))(23) 



and for the normal component: 

[uj + iO+ - — - ^jUqiks - q, ki, tj) = i(5k,ki 

771 ZTO 

+ — ^"'^ 7'-+'^^ (E ^^q(k-, k„ C) + ^q(k., k, , C) 

+Wq(kg~q, ki,c^))(24) 

The prime in the sum excludes terms involving the con- 
densate mode. This close set of equations is solved in the 
Appendix A. For k / kg, kg — q, the scattering solutions 
are lo = Ck+q — ^k- For excitations involving a super- 
fluid mode k = kg,kg — q, the presence of interaction 
with the macroscopic condensate transforms the scatter- 
ing solutions into collective solutions of the discriminant 
equation: 



A(q,cj) = /C„(q,Lj)[(tJ + iO+ - 



kg.q 

TO 



2 B' 



-(^„(q,.)-l)^S!.0 (25) 
mv TO 



C2q2 + (^)2 

Im 



is the Bogoliubov excitation energy, 



47rankj 



is the sound velocity and 



^«(q, w) = 1 



mV ^ UJ 

k 



k-l-q 



(26) 



(27) 



(28) 



is the dynamic dielectric function of the normal fluid. 
Ea. H25(l allows to find a dispersion relation for any func- 
tion 71k possibly in non thermodynamic equilibrium. In 
this sense, it generalizes the dispersion relation that is 
obtained equivalently from the density fluctuations re- 
sponse formalism for a gas at equilibrium 24,, 31, 37, 3i, 
Isgf . Note that in the non number conserving Beliaev 
formalism has been used to derive the dispersion relation 
up to the next order beyond the Bogoliubov theory. The 
solution can be put in the complex form: to — i^q 
where 7q corresponds to the Landau damping. For the 
particular case of a weakly depleted Bose gas, we can 
solve analytically Ea. H25|l . In that case, we can approxi- 
mate ml: 



/C„(q, w) ~ 1 + iIm/C„(q, u) 



(29) 



Ea. (|25() is obtained considering in a first approximation 
that the imaginary term can be neglected. We find that 
the real part corresponds to the Bogoliubov spectrum: 



ks-q B 

q 



(30) 



The imaginary part corresponds to the Landau damping 
and can be calculated perturbatively assuming |7q| -C Wq 
which is the case for a weakly depleted condensate. We 
find, up to the first order, 



7q ~ Im/C„(q,LJq) 



47ranks q 
mV m 



(31) 



By analogy with plasma physics, we say that a Bose gas is 
stable provided that 7q > and unstable otherwise. In a 
stable condensate the collective oscillations are damped, 
while in a unstable condensate they are amplified expo- 
nentially. 

For the case of thermal equilibrium H58I) with Vn = 0, 
/i = and a temperature close to zero, ImA;;„(q, Wq) is 
calculated in appendix B and is a positive function for 
Wq > and negative otherwise. Thus, from H31f) the 
stability condition could be written as: 



TO 



> 



(32) 
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We can check that this inequahty is fulfilled for all val- 
ues of the momentum at the condition that |ks/m| < c. 
In other words, the condensate is stable if its velocity 
relative to the normal fluid is much less than the sound 
velocity. Anticipating the next subsections, this condi- 
tion corresponds to the Landau criterion for superfluidity 
for a weakly depleted Bose gas. More generally, the oc- 
currence of instability depends on the form of which 
influences the spectrum obtained from H25() . In the case 
of a plasma, this problem has been studied long time ago 

nisi. 

Finally, in the absence of a macroscopic condensate i.e. 
rik^/V — > 0, we recover the two scattering solutions for 
the superfluid: 



CJ — ± 

m 2m 

and the collective solution is given by 
/C„(q, w) = 



(33) 



(34) 



According to [31|, at equihbrium, this collective solution 
contains only an imaginary part and thus no collective os- 
cillation can be observed in the system. Therefore, any 
collective oscillation results specifically from the conden- 
sation which transforms the scattering solution of the 
condensed mode into a collective solution. 



C. Derivation of the kinetic equation 

From the results of the previous section, we are now 
ready to derive a generalized Boltzmann like equation for 
the Bose condensed gas. We define the spatial correlation 
function as the average (pk'.-qPk.q)- Using Ea. (|14l) . we 
derive the following equation for the correlation function: 

,.d , (k'-k).q q2 

ot m m 

= —77 kk(2 - (5k,k" - ^k-q,k") 

k" 

-'^k+q(2 - (5k, k" - <5k+q,k")] (/Ok',-qPk",q) 

+ [nk'(2 - (5k', k" - <5k'+q,k") 

-nk'-q(2 - (5k',k" - ^k'-q,k")] (Pk",-q/Ok,q) (35) 

This function can be decomposed as a non interacting 
part and an interacting part 

(Pk',-q/Ok,q) ("k + l)flk"5k'-k,q 

+nknk'(5q,o(l - '5k,k') - n-k(5k,k'^q,o + gq(k, k') (36) 

where (7q(k, k') represents the correlation function due 
to the interactions. For q ^ or k k', the non inter- 
acting part follows the Wick's decomposition and since 
the system is homogeneous (cj^Ck') = ^k.k'^^k- For q = 
or k = k', however, it reduces to the quadratic average 
population that we choose to be (n^) = n^. 



Would we have used the Wick's decomposition in that 
case then (n^) = 2n^ which corresponds to non zero 
particle number fluctuations (S'^ni^) = and the total 
particle number would display fluctuations as well. In- 
deed, since (nk^^k') = '^k'T-k' for k ^ k', we calculate 
(S^N) ~ X^k'^k- f^*5 presence of condensation, these 
fluctuations are huge i.e. of the same order of the average 
value [iBI • Consequently, in an isolated system where the 
total particle number is conserved, this unphysical situ- 
ation must be excluded. 

Taking the average over each side of Ea. (|12|l . a com- 
parison with H36|l allows to deduce: 

* = II ^ (5q(k, k') - gq(k - q, k') 

q,k' 

+gq(k',k)-.gq(k',k + q)) (37) 

On the other hand, inserting this definition into Ea. l|35|) 
and using (|37|l . we obtain 

Ot m m 

= Qq(k, k') + ^ V [nk(2 - (5k,k" - 4-q.k") 

k" 

-?ik+q(2 - (5k, k" - (5k+q,k")]3q(k",k') 
+ [n-k'(2 - (5k',k" - (5k'+q,k") 
-nk'-q(2 - ,5k'.k" - (5k'-q,k")]5q(k,k") (38) 



where 



Sttc 



Qq(k,k') = :^^[('^k - ?^k+q)(?^k'-q + l)?^k' 

-f (nk' - nk'-q)(nk + l)«k+q] 



mV 



('5k,ks'5: 



k',ks ~ <5k,ks-q'5k',ks+q 



[(1 + ?^kj"k,+q + (1 + nk3-q)'T-ks] (39) 

is the inhomogeneous term. To get H39|l . we eliminate 
terms involving delta functions which will give negli- 
gible contribution to the QKE in the thermodynamic 
limit. Both Ea. (|37() and Ea. H38() form a close set in 
which (7q(k, k') must be eliminated in order to get a ki- 
netic equation for nk. This elimination is done following 
an analog procedure to that of Ichimaru We as- 

sume that the correlations due to the interactions are 
non-existent for t — > — oo. This requirement is usual in 
kinetic theory and allows to provide the following initial 
condition (7q(k, k')|(^_oo = 0. Also, we assume that rik 
evolves on a much more long time scale than the duration 
of a collision (7q(k, k') and so is considered as constant in 
solving Ea. H38|l . This approximation amounts to claim- 
ing that the binary collision process is instantaneous in 
comparison with the time associated to the relaxation 
of the system. Inspired by the previous subsections and 
by 40], we can check that the solution, satisfying both 
Ea. (|38|) and the initial condition, is expressed in terms of 
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the dielectric propagator as: 

/"OO 

g^{k,k\t) = -i / dt' 

Wq(k, kl , t')W-q(k', k'l, i')Qq(kl, k'l , < - i') (40) 

We have inserted the exphcit time dependence in the 
functions. If the creation of such correlations is much 
faster in comparison with the relaxation time for n^, 
(5q(ki, k'l, t — i') ~ (5q(ki, k'l, and we obtain a Marko- 
vian equation. We re express the dielectric propagator in 
Fourier transform according to, 

iYq(k,ki,i) = / dwe-'"%(k,ki,c^) (41) 

where the contour C extends from —oo to +00 along a 
path in the upper half of the lu plane in such a way that 
all the singularities lie below it. The substitution allows 
to carry out successively integrations over t' and uj' by 



closing the contour in the upper half plane in order to 
get: 



3q(k,k')- 




Wq(k, kl, C^)iY_q(k', k'l , -L0)Q^{kl , k[,t) (42) 

Calculations in appendix C allow to find an explicit ex- 
pression for gq(k, k') in terms of the one particle distri- 
bution function, provided the Landau damping factor is 
always positive. The substitution into l|37|) allows finally 
to get the GRPA kinetic equation for a stable Bose gas: 

d 

= [nk';ks] = Ck[nk';ks] -l-Ck[nk';ks] (43) 

where we define the collision terms as a functional of rik' 
and a function of kg . The first term describes the collision 
rate between particles of the normal fluid: 



Ckb^k'; kg] = ^ 

q,k' 



Sira 
mV 



^(q, Ck+q - Ck) 



(1 — ^k.K ^ '5k+q,ks ^ <5k',ks — <5k'-q,ks)7''^(ek+q + Ck'-q — Ek — Ck' ) 

[nk+q«k'-q("-k + l)(?lk' + 1) - "-k"-k' ("-k+q + l)("-k'-q + 1)] (44) 



The second term describes the collision rate between condensed and non condensed particles, the condensed particle 
is cither the input or the output state in the scattering process: 



Ck[nk';ks 



E 



\mV} 



('5k,ks + (^k+q.k, + Sk'M^ + <5k'-q,kJ 



q.k' |^*(q, Ek+q - ek)^(q, fik+q - ^k) 
Tr6{ek+q + Ek'-q - Ck - ek')["-k+qflk'-q("-k + l)(?^k' + 1) - "-k"-k' ("-k+q + l)(nk'-q + 1)] (45) 

I 



In the first expression H44|l. the contact potential has been 
replaced by an effective one depending on the transfer 
particle energy u — Ck+q — e^: 



Sira 

IC{q,uj) 



(46) 



The correcting term is the dynamic dielectric constant: 
/C(q,tj) 



7nV m 



(uj + iO+ - - (^31)2 + l^TOTM^q^ 



A(q,w) 



(47) 



(^ + iO _ k^)2 _ (^)2 + 4^q2 

\ ' ^ 7n ' ^ 2m -' mV m 

In the second expression H45(l . we define another dynam- 
ical dielectric function: 



/C(q,cj) = 



A(q,c^) 



(48) 



This result generalizes in the case where ^ 0. 
Plugging this expression into Ea. H45|l . the energy con- 

k ^ 

servation imposes the two choices uj — ± when 
a particle of the condensate participates to the collision 
process. But, with such an energy transfer and for a 
macroscopic population of the condensate, this dynamic 
dielectric function gets infinite: 



^(q'^)l,.,-iss 



(49) 



leading to Cki^^k'! kg] — > 0. In this situation, surprisingly, 
no collision involving the condensate particle occurs. As 
a consequence, in a uniform Bose gas, no transfer of par- 
ticle is possible between the condensate and the normal 
component. Indeed, due to the absence of Fock inter- 
action term in the GRPA, the effective potential has a 
different expression when a particle of the condensate is 
involved in a scattering. As said in section 2, a potential 
induced by the condensate compensates exactly the po- 
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tential created by the excited particle susceptible to scat- 
ter. As a consequence, the dielectric function suppresses 
the effectiveness of the potential which thus becomes in- 
visible to the condensate and totally shelters it from col- 
lision. Furthermore, let us note that this shielding re- 
mains unchanged for any momentum kg which preserves 
the stability of the condensate. If the Landau damping 
factor 7q becomes negative, then an instability occurs 
and therefore Ea. H44|l and Ea. H45() are not longer valid, 
as well as all considerations concerning collision block- 
ade. In that case, a more elaborated derivation must 
be carried out that could allow particle exchange with 
the condensate. The expressions HC6|I and (|C7|I must be 
recalculated taking into account the instability [sS, ■ 

In this way, the GRPA kinetic equation provides a 
different understanding of the superfluidity phenomenon 
as it describes the motion of two independent fluids 
that cannot transfer particle through collisions. Since 
^rikj, = 0, is really a constant of motion indepen- 
dent of the dynamics of the normal fluid. 

The collision blockade is the consequence of a higher 
order expansion in the interaction parameter. Since 
/C(q, a;)|a^o — 1, up to the second order in the inter- 
action, we recover the UUQKE allowing collisions with 
the condensate. Only an infinite re-summation of ap- 
propriate higher order contributions through GRPA gen- 
erates collectively a perfect collision blockade. Also, in 
the limit where the condensate is not macroscopically 
populated ^ks/^ — > 0, then IC{q^,uj) /C„(q, w) and 
/C(q, A^„(q, a;) and collisions with particle of the 

mode kg become again possible. The Ea. (|^ is a dif- 
ferent version of the kinetic equations formulated by 
Zaremba et al. 2] in which collisions are possible with 
the condensate particle. Indeed, in their approach, the 
authors have done an approximation which is equivalent 
to replacing the energy ek by IjKil) and to suppressing the 
integral terms of the right hand side in (|14l) . As a con- 
sequence, they obtained a QKE in which the dielectric 
function does not appear and with the feature of having 
a gap in the transition of a particle to the condensate. 

Let us notice that the same phenomenon occurs for the 
simple RPA where the exchange terms have been omit- 
ted. In the SRPA the calculation is similar to those lead- 
ing to the quantum kinetic equation for a plasma. For 
details of such derivation, see Balescu j^^. The only dif- 
ference is that the coulombian potential becomes a con- 
tact potential and that the fermion now becomes a boson, 
meaning that any factor 1 — rik now becomes l-l-rik in the 
collision term. In that case, the dielectric function has 
the same form for both the condensed and the normal 
modes Isil: 



IC 



SRPA 



"k - nic+q 



k ^ 

47rani, 



iOj 



k.q 



31 

2m 



mV 



The normal dielectric component A^f^^^(q, w) is of the 
same form as Ea. (|28f) . except that the exchange inter- 
action term has been removed. Similarly, the particles 
of the condensate are protected from any collision with 
the surrounding. Nevertheless, in the SRPA, collisions in- 
volving low level excited states are strongly prevented be- 
cause the energy transfer uj is close to the one correspond- 
ing to the condensate. On the contrary, in the GRPA, 
a supplementary term coming from the exchange effect 

ensures a reasonable value for the effective po- 
rn v ^ 

tential guaranteeing an efficient collision rate, even with 
particle in the energy levels close to the condensate level. 



D. Properties of the RPA collision term 

The collision term exhibits a number of remarkable 
properties analog to those encountered by other kinetic 
equations j4ii] . These allow to establish the particle num- 
ber, momentum and energy conservation laws as well as 
the Boltzmann H-theorem. Since the condensed particles 
do not participate to the collision process, these concerns 
only the excited particles of the normal fluid. Three of 
them can be stated as: 



^Ck[nk';ks] = 

k 

^kCk[7ik';ks] = 

k 

I 

^ekCk[nk';ks] = 



(51) 
(52) 
(53) 



These three integral relations are checked easily by di- 
viding the collision term in four equal parts. After 
carrying out the successive changes of the integration 
variables k ^ k', q ^ — q on the second term, k <-> 



q, k' 



k -I- q on the third term and a succes- 



sive combination of these two variable changes in the 
fourth term, these four terms cancel each other if we 
use the relation /C(q, w) = 1C*{—ci^—lo). As a conse- 
quence, for a uniform Bose gas, the total particle number, 
the total momentum and the total kinetic energy, asso- 
ciated with the normal component are independent of 
the time i.e. [d/ dt)^^n'^ = 0, [d/ dt)^^'kn'^ — and 
(d/dt) ^jj. ek?^k — 0. These properties together with the 
conservation of n\^^ imply that the total energy in the 
Hartree-Fock approximation is also conserved [3 i.e.: 



dt^^ 



ek"-k 







(54) 



Another crucial property is the Boltzmann H-theorem. 
The entropy, due to thermal excitations, has the expres- 
sion: 



(50) 



S 



E 

k 



[(1 



lln(l 



. In n(. 



(55) 
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The time evolution of the entropy is always positive. A calculate a positive production of entropy: 
similar derivation using the change of variables allows to 



dS_ 
It 



= E 

q,k,k' 



Sira 
mV 



^(q, Ck+q - Ck) 



-SiCk+q + Ek'-q - Ck - Ek') In 



<+q"k'-qK + l)K'+l) 
"k"L'K+a + l)K'~a + l) 



4'-qK + 1)K' + 1) - «'K+q + l)K'-q + 1)] > 
I 



(56) 



The positivity is a consequence of the mathematical re- 
lation \ii(x/y){x — y) > 0. The equality is achieved for 
X = y. From Ea. H56|) we deduce that the entropy always 
increases until the system reaches a stationary equilib- 
rium distribution. This occurs when the production of 
entropy becomes zero. In that situation, 



k+q"^k'-q 



'k'^k' 



K+a + l)K'-q + l) K + l)« + l) 



(57) 



This relation holds only for the Bose-Einstein distribu- 
tion function jlj|: 



(58) 



exp [/3(ek - k.v„ - ^)] - 1 



where Vn is the average velocity of the normal compo- 
nent. In this way, the inverse temperature [3 and the 
chemical potential /i arc defined as the free parameters 
of the equilibrium solution. Although similar, the prop- 
erties of the thermodynamic equilibrium predicted by the 
H-theorem in the RPA are really different from the one 
predicted by the calculation of the ensemble partition 
functions of an ideal Bose gas The superfluid and 

normal components can move relatively with two differ- 
ent velocities Vg and Vn. The relative difference Vg — Vn 
and /i are subject to the only constraint of stability which 
gives a limitation. As said in the previous subsection, 
for Vn = 0, and the temperature close to zero, and us- 
ing (jB4|l in (|31|l , this constraint is realized provided that 
jjL — Q and that the relative velocity does not exceed the 
sound velocity (see Eq.|(211)). Otherwise, the condensate 
becomes instable. Thus, we recover the Landau criterion 
for the weakly interacting Bose gas. Note that this result 
differs from the ZNG approach since it corresponds to an 
equilibrium fugacity of the normal fluid equal to one (see 
Eq.(44) in :2|). 

On the contrary, the statistical equilibrium ensemble 
formalism for an ideal Bose gas imposes a zero relative ve- 
locity and a chemical potential close to zero when conden- 
sation occurs 15]. This contradiction might be explained 
if we remind that in this formalism we postulate an equal 
a priori probability of any possible configuration of the 
gas |l2l [l5l |. This basic assumption of equilibrium sta- 
tistical mechanics originates from the observation that, 
over a long time range, the collision process will mix sta- 
tistically these configurations. But, since in the GRPA 



kinetic equation the collision is blocked, the equal a pri- 
ori probability might not work as far as the condensate 
mode is concerned. This observation suggests that the 
condensed particle population for an homogeneous sys- 
tems is not distributed anymore but has a well defined 
value. 



E. Condensate formation 



The kinetic equation (|43|) does not explain the conden- 
sate formation through, for example, evaporative cooling 
[33.l35ll3^ . If we start indeed from an initially condensed 
gas with a Bose-Einstein momentum distribution whose 
the tail has been cut, then according to the RPA model, 
the irreversible evolution process towards equilibrium will 
undergo a strange behavior. Instead of populating the 
lowest energy level which is forbidden, the excited par- 
ticles will be transferred to the lowest excited levels in 
which a macroscopic population can eventually appear. 
In that situation, we are led to a fragmentation of the 
condensate into at least two energy levels. But such two 
macroscopic states contradict the assumption made at 
the beginning that only one macroscopic state does ex- 
ist. A more elaborated model can take into account more 
than one level macroscopically populated. Although this 
more general description must not be excluded, for ener- 
getic reason, this phenomenon is not likely to happen . 
A fragmented condensate has a much higher potential en- 
ergy than a non-fragmented one due to the presence of 
a Fock interaction energy between two fragmented parts. 
To explain the condensate formation and if we exclude 
fragmentation, we must assume either that the Bose gas 
must be instable due to a far non-equilibrium situation or 
that it must be inhomogeneous with regions of strong de- 
pletion. Indeed, in those regions - like the edge of the gas 
- the condensate population is not macroscopic anymore 
and so collisions between the two fluids may happen. 
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IV. EXTENSION TO A WEAKLY 
INHOMOGENEOUS BOSE GAS 

A. Inhomogeneous equations 

The extension of the previous equations to inhomo- 
geneity is important in order to understand how the evo- 
lution of the condensate and the thermal excitations are 
coupled through mean field forces. In what follows, we 
assume that the gas is weakly inhomogeneous i.e. most of 
the quasi-particles collide inside a sufficiently small vol- 
ume that can be considered as homogeneous. To quan- 
tify the level of acceptable inhomogeneity, we divide the 
volume V of the gas into small cubic volume Q. The 
edge I of each volume = ^■^ is adjusted in such a way 
to be infinitesimal so that the gas is homogeneous in- 
side it, but big enough to contain a large amount of 
particle whose dynamic obeys still locally the homoge- 
neous equation. In the literature I is referred as the 
hydrodynamic scale and is estimated from the formula 
l{r,t) = /9e(r, t)/VrPe(r, where pe{j^,t) is the normal 
fluid density and r is the position in space of the small 
volume. This length scale must be much greater than 
the mean free path which can be estimated more less as 
l/(pef) where a = Stto^ is the total cross section 12]. 

In this infinitesimal volume, we define the local con- 
densate wave function 



(59) 



which corresponds to the eigenfunction of the density ma- 
trix with the highest macroscopic eigenvalue The 
local momentum of the condensate depends also on the 
position and the time through the relation 



ks(r,i) ^ Vr6'(r,i) 



(60) 



Up to a constant phase, the amplitude and the local mo- 
mentum of the condensate characterize fully '^(r,t). In 
this volume, we define also the local particle number dis- 
tribution or Wigner function njj(r, t) (for clarity, explicit 
dependence of time has been added). Also the particles 
feel an external local potential Ve//(r). Up to the first 
order in the interaction potential, the local energy den- 
sity is divided into the kinetic part, an external potential 
part and the Hartree-Fock part |3, 0| '■ 

+ E(£ + ^e//(r))<(r,t)-f^[|*(r,t)r-f 

4|vl/(r, i) P E + 2(E (61) 



This energy density varies slowly in space. The gradients 
produced by these variations govern the dynamic of the 
gas particle between each small volume ft. The local 
energy per excited particle with a momentum k is given 
by the derivative of the local energy density with respect 
to the particle number 0: 



ek(r,0 



(62) 



where we define the effective potentials felt by the con- 
densate and the particle 



K(r,t) 



Vextir) + — 
mV 



l*(r,t)P+ V<(r,i) 



k 



(63) 



Between each infinitesimal volume, particles are trans- 
ferred due to the gradients of the particle density and 
of the potential energy. As a consequence, the kinetic 
equations must be modified in order to take into account 
locally the modifications inside O. If we consider that 
the excited particle moves classically, then the dynamic 
of transfer is given by the Liouville operator acting on 
the distribution function: 



{ek(r,0,<(r,i)} 
= Vkek(r, i).Vk<(r, t) - Vrek(r, t).Vk<(r, t) (64) 

The condensate particles however move locally accord- 
ing to the quantum wave function, which obeys to the 
generalized Gross-Pitaevskii equation with the effective 
potential j^j: 



V^c(r, t) 



, . Ana 
^-(r) + ;^ 



\^ir,t)\' 



2^<(r,t) 

k 



(65) 



Furthermore, we assume that locally in the small vol- 
ume fl: 1) the system is sufficiently homogeneous so that 
the expressions (|44ll and H45|l remain valid and the con- 
densate is not affected by collision; 2) the condensate is 
locally macroscopically populated. Combining these re- 
sults together with the collision terms, we can write the 
kinetic equations for weakly inhomogeneous and stable 
Bose gas in the region of condensation. We find two cou- 
pled set of equations, one for the condensate, the other 
for the normal fluid: 
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9, 



2m 



+ K(r,t) 



*(r,t) 



dt 



nk(r,t) = 



.Vr + VrV;(r,i).Vk 

ni 



4(r, i)+Ck[nk'(r,t);ks(r,<)] 



(66) 
(67) 



With the exception of the colhsion term, the Eg. 166(1 
and H67|l are identical to the kinetic equations formu- 
lated by Zaremba et al. The inhomogeneous ki- 
netic equations satisfy locally the conservation laws in 
these regions of highly populated condensate. We can 
derive indeed a conservation equation for the local par- 
ticle number N{r,t) — X]k"'k(r,i), the local momentum 
P{r,t) = ^,^knk(r,t) and the total local energy £{r,t) 
It is also easy to verify that the pro duction of the 
local entropy S{r,t) is always positive [ia |. 

The local production of entropy stops when we reach 
the local equilibrium: 



feq 



(r,t) 



exp [/?(r, t){ek - k.v„(r, t) - ^(r, t))] - 1 



(68) 



where /3(r, i), Vn(r,t) and ^(r,i) are now local functions 
of the position and the time. 

A gap shows up between the energies of the condensed 
and non condensed particles since Vc{r,t) < Ve{r,t). 
This is not a problem as long as in the region of the 
gap the transfer of particle is forbidden. The transfer of 
particle is only possible in the region where there is no 
gap i.e. when the condensate population is not macro- 
scopic. Therefore, on the basis of these non-equilibrium 
considerations, the Hartree-Fock model showing up this 
forbidden gap does not suffer from any inconsistency re- 
lated to the conservation of energy. 



B. The superfluid universe at finite temperature 

Using the expression (|68|) as the solution of the ki- 
netic equation, we can find a set of equations describing 
the non dissipative motion of the condensed gas at fi- 
nite temperature. If we assume that Vn(r,i) — then 
the substitution of H68|) in Ea. (|67|l imposes that n'^'^{r, t) 
must be stationary in time, that /3(r, i) = /3 is a constant 
and that fJ.(j,t) = —Ve{r,t) + /ig, where He is a chemi- 
cal potential independent of the position controlling the 
number of excited particles. As a consequence, after car- 
rying out the integration over the momentum, the local 
number of excited particle is given by the closed equation: 



Ar|'?(r)=^<«(r) 



„/3[M.-Vext(r)-|jii(|*(r,t)|^+JV|''(r))] 



where gi{x) = J2'jLi j^^'' ■ Since the system is station- 
ary, the macroscopic condensate wave function has the 
form: 



*(r,i) = e-'^'^*^'c(r) 



(70) 



where fj,c is the chemical potential of controlling the pop- 
ulation of the condensate particle. The substitution of 
this form into Ea. (|66|l produces: 



V2 
^ r 

2m 



A.'jrn 

+ K.t(r) + — (|*,(r)p + 2N!'>{r)) 



mV 



*c(r) 



Aic«'c(r)(71) 



The coupled set of equations and ifTT)) describe lo- 
cally all the non dissipative or superfluid phenomena. 
These are valid provided the gas is stable. As an ex- 
ample, they describe the superfluid moving with a dif- 
ferent velocity than the normal fluid. In that case when 
Ve2:t(r) = 0, the solution is a plane wave function: 



*c(r) 



(72) 



(69) 



and N^'>(y) is a constant. More generally, any non dissi- 
pative complex structure at finite temperature, like vor- 
tices, should be a solution of H69() and H71|l . 

In contrast to previous studies, non-dissipative phe- 
nomena are not a requirement or assumption of a model 
but rather a prediction of the GRPA kinetic theory. 
Another difference is that the chemical potential for 
both fluids /ic and /ig must not be necessarily iden- 
tical, contrary to the equilibrium statistical mechanics 
which imposes equality Jl5l|. In our example, = 
{^T:aN)/{mV) -/-[ic = ks7(2m) -h47ra(2iVe-hnkJ/(my) 
at the difference of where [i^ = /^c for kg = 0. 



V. ANALOGY WITH PLASMON THEORY 

The derivation leading to the QKE for a Bose gas is 
similar to the one leading to the QKE equation for plasma 
physics 'Q&^ J£, 40J except that, in these equations, the 
collective excitations are derived in the SRPA. In plasma 
physics, the Coulomb interactions potential for the col- 
lision process is screened beyond the Debye wavelength. 
The dynamic dielectric function displays this screening 
and removes the singularity in the Coulomb potential in 
the long wavelength limit. Static screening is well known 
since the Debye theory. However, for dynamic screening, 
Wyld and Pines have proposed an interpretation in terms 
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of the plasmon theory. This theory is an alternative ver- 
sion of the kinetic theory of a plasma which emphasizes 
the role played by the collective modes. According to 
their work, the effective potential interaction is the result 
of a plasmon mediating the interaction with the quasi 
particles. In other words, during a collision process, a 
quasi-particle emits an intermediate plasmon on the en- 
ergy shell with the momentum transfer q. Later on, this 
plasmon is eventually absorded by another quasiparticle 
which acquires a new momentum and a new energy. The 
energy spectrum for this plasmon corresponds precisely 
to the frequency spectrum of the collective excitations in 
a plasma. 

In the case of a Bose condensed gas, the collective ex- 
citations spectrum corresponds to the Bogoliubov fre- 
quency spectrum for low temperature. This property 
suggests the interpretation that these phonon-like exci- 
tations play also the role of mediators during collisions 
between quasi-particles. Thus, the Bose condensation 
has the effect to transform the scattering modes, used 
for collision involving condensed particle, into a collec- 
tive mode used for mediating the interaction between non 
condensed particles. 

The plasmon theory has been derived using the theory 



of quantum electrodynamics. In principle, for a Bose gas, 
a similar approach must be carried out taking into ac- 
count that, at a more fundamental level, the interaction 
potential originates from processes involving the absorp- 
tion and the emission of photons. In this paper, we shall 
not rederive a similar theory but rather recover heuristi- 
cally the analogy with the plasmon theory. 

Following this approach, this intermediate excitation 
behaves like a particle characterized by its own distribu- 
tion function /q and that we shall call by analogy "con- 
denson" . The energy spectrum of the condenson is given 
by the zeroes of the dynamic dielectric function i.e. by 
the solution of A(q, w) = 0. We consider the simple case 
of an homogenous condensate at rest kg = and weakly 
depleted. In that case, the solution is the complex num- 
ber to = LUq — Z7q, whcrc the real part represents the 
energy spectrum tjq — and the imaginary part repre- 
sents the decay rate of the condenson given by Ea. ljHlll . 
Then, looking at Eq.(7) and Eq.(8) in .20^, we can write 
two coupled equations, one for the dynamic evolution of 
the quasi-particle and the other for the evolution of the 
condenson. They describe the time rate change of these 
quantities due to the emission and absorption of a con- 
denson by a quasi-particle: 



q ^ 



+ ^<5K - ek + ek-q) (« + iX-q/q - <«-q + l)(/q + 1)) ] (73) 

9 r r Stto 87ranoq^ ^, x / / / -, x 

-/q = -27q/q+— ^-^^^5K-ek + ek-qKK-q+l 74 

at mV m'^VuJa -i >d 



Similarly to the reasoning of Wyld and Pines, the ki- 
netic equation (|44|l can be rederived from Ea. H73(l and 
Ea. (|74|l by eliminating adiabatically /q. We assume that 
the time derivative in Eq. (|74|l is zero and the elimination 
allows to obtain an equation for n'^ only. Ea. l44() is re- 
covered, provided that the condensate is weakly depleted. 
Under these circonstances, the square of the effective po- 
tential has a narrow Lorentzian shape whose two peaks 
and widths correspond to the condenson energy and de- 
cay rate respectively (see (22j]). It plays the role of a 
propagator for the mediated interaction. This supposes 
that the dielectric function can be approximated around 
the resonant frequencies as (see Eq.(25) in [2^ for com- 
parison): 



leading to the approximation: 



1 



2\ 2 

{5{lO - Wq) + 5{U0 + Wq)076) 



TT / 47ranoq' 
|/C(q,t^)|2 7q \ m^Fwq 

These coupled equations suggest that a weakly interact- 
ing Bose condensed gas is in reality composed of two kind 
of bosonic particles: the quasi-particles that become the 
real particles in the limit of an ideal Bose gas and the 
condensons that result from a mediation (see Fig.l). At 
thermal equilibrium, the various particles obey to the 
Bose-Einstein statistics: 



1 



exp [/3(ek - m)] - 1 



(77) 



and 



/C(q,tj) 



47ranoq-' 
47ranoq"^ 



.(c 



^q 



«7q), 

- i7-q) 



'^q 



/q = 



exp (/3a;q 



1 



(78) 



(75) III this way, the condenson distribution function has the 
1 same form as the excitations distribution function pre- 
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FIG. 1: Feymann diagram of the interaction of two quasi- 
particle (full line) mediated by the condenson (wavy line) 



dieted by Bogohubov 18]. But there is an important 
difference however: 

In the BogoHubov theory, the Bogoliubov excitations 
correspond to the quasi-particle. On the contrary, in the 
GRPA model, the Bogoliubov excitations correspond to 
the condensons and, in this respect, are not the quasi- 
particles. 

The QKEs proposed by 0, H, 0, US suggest instead 
that the Bogoliubov excitations correspond to the quasi- 
particle. We must be cautious with this interpretation 
since the GRPA model is not an exact one. It might be 
that, to a next order expansion in the interaction, the 
quasi-particles have a different energy spectrum that is 
hnear in the momentum. Without any further investi- 
gation, we cannot claim that the Goldstone boson, that 
would be normally predicted from the U{1) theory, is the 
condenson or the quasi-particle. 

The condenson has a decay rate 27q. Thus, a duration 
of the collision is limited by its lifetime and can be con- 
sidered as instantaneous, if much lower than the inverse 
relaxation rate. Assuming that the particle has a mean 
velocity l/(/3m)^/^, this requirement imposes for the typ- 
ical values: 7q ^ pecr/{(3my^^ . If it is not fulfilled, then 
the long lifetime of the condenson suggests that its dy- 
namic evolution is not instantaneous and that the full set 
of equations (|75|) and l(7I|) must be rather used. 



VI. EQUIVALENCE WITH THE BOGOLIUBOV 
THEORY 



The RPA model presents important differences com- 
pared to any previous model for finite temperature. How- 
ever for zero temperature, it allows to recover the next or- 
der correction for the ground state energy and the static 
structure factor usually obtained from the Bogoliubov 
theory. Such a derivation has already been made by 
Nozieres and Pines for an electron gas using the den- 
sity response formalism. They were also able to apply 
this formalism to the condensate The nice result of 
this section is that, precisely, the formalism developed in 
this paper allows to recover exactly the same results for 



the ground state energy and the structure factor. 

Indeed, the correlation function defined in H36|) allows 
to calculate the correction to the energy but also to the 
total particle momentum distribution N^. This new dis- 
tribution differs from the quasi-particles population dis- 
tribution nk. We start initially from the distribution 
without interaction = Sk oUo where only the con- 
densed mode kg = is populated. Then, we switch on 
adiabatically the interaction and the correlations appear 
in a form given by (^U)). As a consequence, the conden- 
sate becomes weakly depleted and the total kinetic and 
potential energy gets modified. 

These corrections can be calculated directly the static 
structure factor which, using and for q 7^ 0, is given 
by: 

NSiq) = (pIp^) = ^(nu + l)nk+q + J2 5q(k, k')(79) 

k k,k' 

where N — X^k ~ Sk'^k- An expression for this 
integral can be found in the appendix C. In contrast, 
Nozieres and Pines express the structure factor in terms 
of the susceptibility function © of the ground state as 



NS{q) ^ (ptpq) 



div 



Imx(q,w) 



(80) 



At zero temperature, the limit ri[^ must be taken 
carefully. It implies that the Landau damping 7q ^ 
but this parameter is different from 0+ . The limit 0+ — > 
must be carried out before the limit 7q — *■ and this pro- 
cedure can be done by adjusting adequately the infinites- 
imal value of 0_|_ and n'-^. In particular, for a normal fluid 
in equilibrium, this amounts to prescribing to take the 
limit 0+ to zero and then after taking the limit /3 — > 00 
in Ea. (|58|l . The physical reason for such a procedure is 
that the presence of infinitesimal fraction of excited par- 
ticles will force the system to create the correct equilib- 
rium correlation. In these conditions and using Ea. H29|) 
and Ea. H31() . we can re express for kg = 0: 



(81) 



Also, we note the relation: 



lmA'{q,uj) 



mV 



-Im/C„(q,c^) (82) 



Sira exp{(3uj) — 1 

Inserting these results and ll^ into (|C13I) . setting to 
zero the real part of /C(q, w) and A'{q,uj), and setting 
Ao{q,uj) = N/[LJ + iO+ + q^/{2m)] then Eq.lO becomes: 

2 



^(q) 



1 - 



duj' 



\2m J 



iTraN q 
mV m 



q2 g/3w _j_ ^ ' 



7q 



-> 
00 



(83) 
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In order to get H83|l . the integration has been carried out 
transforming the Lorentzian factors into delta functions 
and using (jSIJ. This result is identical to that obtained 
from the Bogoliubov theory. This result can also been 
obtained substituting ^ into H80() and carry out the in- 
tegration over uj. Instead of evaluating the total energy 
and the correction to the particle number distribution 
directly, some useful tricks allow to avoid some complex 
calculation . The total ground state energy is a func- 
tional of the interaction a and the dispersion relation 
Ek (in what follows, we consider only the first term in 
It can be expressed in terms of the matrix element 
i?(a,ek) — {^P\H\iIj) where \tp) is the ground state de- 
pending also on a and ek- Using the property that 
is a normalized eigenfunction, the first derivative with 
respect to these parameters produces: 



a = £'int(a,ek) 

9£:(a,ek) 



dek 



(84) 
(85) 



The total interaction energy can be expressed in terms 
of the structure factor: 



get for k 7^ 0: 



iVk 



and for k = 0: 

No = N - 



1 / i£_ J- 4;rajV 
I 2m ^ mV 



(90) 



(91) 



From our number conserving approach, we recover the 
well-known results usually obtained from the Bogoliubov 
approach at zero temperature, i.e. the first order energy 
correction in H89|l and the first order correction to the 
number of excited particles (|90|l showing a depletion of 
the condensate particle. In principle, the present method 
can be extended to finite temperature, since the relations 
and (|SS|l are valid also for any excited state labeled 
by the occupation number rik. The only problem remains 
the explicit calculation of (PqPq)- 



VII. CONCLUSIONS AND PERSPECTIVES 



27ra 

Eint{a, Ck) = — 77 
mV 



7V(7V-l) + iV^(5(q) 

q#0 



(86) 



At zero temperature and without interaction, all the par- 
ticles are in the ground state and thus E{0,e\^) = A^eo- 
With this initial condition, an integration over a allows 
to calculate the total energy: 



E{a,ek) = iVeo 



° da' 



Eint (a', Ek) 



V 



(87) 



where a supplementary term coming from ^ has been 
added to remove ultra-violet divergencies. The structure 
factor can be recalculated for any arbitrary dispersion 
relation ek- Plugging this new expression into H86|) and 
carrying out the integral over a, we obtain: 



£;(a,ek) =iVeo 



27raiV2 
mV 



2 ^ 

q#0 



87ra7V ^ 



,27raA,T m 



A-KaN 



mV 



where Aca 



(eq + 



2eo)/2. In the usual case. 



Agq — q^/2m and an integration over the momentum q 
allows to find the ground state energy [l5| : 



E{a) 



2TTaN^ 128 
mV ^ Is" 



(89) 



Finally, the derivative with respect to ek gives the first 
correction to the momentum particle distribution. We 



We attempted to reexamine the kinetic theory of the 
weakly interacting and stable Bose condensed gas with 
the objective to explain superfluidity for any tempera- 
ture and, ultimately, to be able to distinguish between its 
dissipative and superfluid behaviors. A different QKE, 
taking into account some higher order terms in the inter- 
action parameter, has been derived from the microscopic 
theory and has the merit to predict a collision blockade 
between condensed and non condensed particle. More 
precisely, the condensate remains invisible to any ther- 
mal quasi-particle due to an induced force, which acts as 
a " coarse- graining ", removing any local force necessary 
for binary scattering. In this way a metastable state can 
be built locally with a non-zero relative velocity. For 
the homogeneous gas, a microscopic derivation has been 
carried out using the same expansion procedure as that 
for deriving the QKE for a quantum plasma. The only 
difference is that our approach takes into account the 
exchange term which is of the same order of magnitude 
than the direct term. For the weakly inhomogeneous 
gas, the derivation needs some supplementary justifica- 
tions from first principles, in particular, for the way to 
obtain the generalized Gross-Pitaevskii equation from a 
number conserving approach, but also for the limit of use 
of the collision terms H44|l and (|45|l . 

The generalized RPA is the basic approximation which 
is at the origin of the prediction of the collision blockade 
phenomenon. It is not only a method to achieve results 
but has also the deep meaning that average contribu- 
tions that are not oscillating in phase can be neglected 
for diluted gas. The RPA does not serve only to predict 
the QKE; it allows to build a genuine alternative the- 
ory that can be compared with all previous approaches. 
The main advantages of this theory are: 1) it is number 
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conserving; 2) the statistical equilibrium is not imposed 
as a postulate but is deduced from the QKE and allows 
superfluidity; 3) it is valid for any range of temperature 
below and above the transition point. 

The predictions calculated from this theory, namely 
the collective excitations, could be compared with re- 
sults predicted mainly by the Popov theory or by the 
Hartree-Fock-Bogoliubov theory in the Popov approxi- 
mation [20I . The RPA theory allows also to deter- 
mine the correlation function (;q(k, k') with which one 
can calculate the next order correction. In particular, 
we recover the Bogoliubov results for the next correction 
to the ground state energy and for the particle momen- 
tum distribution. In principle, these corrections can also 
be determined at finite temperature and compared with 
other previous works. Another interesting quantity is 
the static structure factor S{q) which can be evaluated 
directly from the correlation function. 

The RPA theory for dilute Bose gas has some resem- 
blance with the plasma kinetic theory. Both predict that 
scattering proceeds by means of an intermediate excita- 
tion. For a plasma, this excitation is a plasmon responsi- 
ble for the collective plasma oscillation. For a Bose gas, 
it is a phonon-like excitation with a Bogoliubov energy 
spectrum at zero temperature and that, by analogy, we 
can call "condenson" . 

Let us also mention the limitation of the present the- 
ory. Firstly, this theory preserves the existence of a gap 
in apparent contradiction with the Hugenholtz-Pines the- 
orem J8J. We consider only the Hartree-Fock approxi- 
mation for the potential energy of the quasi-particle, ex- 
cluding in this way higher order correlated term. This 
contradiction might be solved by claiming that the poles 



of the one particle Green function have a richer structure; 
namely one gapless pole which represents the Bogoliubov 
or condenson excitation spectrum and another pole which 
may have a gap and which represents the quasi-particle 
spectrum. In this way, we recover the compatibility with 
the HP theorem. Secondly, in principle, the calculation 
must start with the real potential and not the effective 
contact potential (|2Jl. A resummation of all ladder dia- 
grams should allow to reexpress the effective interaction 
in terms of the T matrix for the binary collision. Then 
the limit of low energy allows to deduce the expression ^ 
Q. If rather the T many body matrix is considered, then 
the effective interaction could depend also non linearly on 
the quasi-particle number. Thirdly, the assumption that 
the gas is weakly homogeneous might not be true, due to 
a sharp trap potential realized in real experiment. The 
theory must be improved by decomposing the one parti- 
cle Wigner function in terms of its eigenfunctions and by 
deriving an equation of motion for each of them |T^ . 

Finally, all the predictions, that the RPA theory might 
suggest, must be ultimately confronted to experiments 
in order to be validated. In particular, the analysis of 
the various modes for the collective excitations and its 
damping must be recovered 42]. 

APPENDIX A: SOLUTION OF THE EQUATION 
OF MOTION FOR THE EXCITATIONS 

The system and can be solved exactly. 

We isolate the dielectric propagator on the left hand side 
of each equation. After summing over the superfluid and 
normal modes separately and using H28() , we obtain after 
rearranging: 
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^Z^q(k',kl,C.) = 



ks.q q% „ , ... kg.q q^ 
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(/C„(q,w) - 1) (Wq(ks,ki,w) -^Wq(ks -q,ki,w)) +X:„(q,w)^Z^q(k',ki,u;) i 



^ m 2m 



(A2) 



where we define the function Ilk = i^k^.k + (^k^-q.k- Using the definition H25|l . the solution of this coupled set of 
equations gives: 



Wq(ks,ki,tj) +Wq(ks - q, ki,Lj) = i- 



/C„(q,..)nk, \{u + ^0+ - ^ f - i^f] + ^^(1 - Hk, 



VM; / L ' ^ rn 2m 



(A3) 



and 



' , (1-/C„(q,c.)) {cJ + ^0^-^^)^^{£.f nk, + [(c. + zO+-%^)2-,f](l-nkJ 
2^Wq(k',ki,w) = i ^ J ^ (A4) 

k' 



18 



Consequently, the definitions (|48|l and (|49() allow to write 

^Wq(k',kl,C^ 



(A5) 



Plugging these results into the r.h.s. of H22|l (|23|l and (|24|l . the integral terms do not appear anymore and we get the 
solution: 



^Yq(k, ki,w) = 



LU + iO+ — 



k.q 

m 2m 



<5k.ki + 
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where we define the reverse dielectric function 

(i-nk)(i-nkj , (i-nk)nk, +nk(i-nkj + nknk, (i-/c„(q,c^)/2) 



JC{(l,uj) 



/C(q,a;) 



(A6) 



(A7) 



This function can be rewritten in matrix notation, where we distinguish in the column and the row the channel kg 
and kg — q from the other channels: 



^k,L(q,^) 



/c-Hq,^) 



/C-l(q,c.) 

la 

2 



(A8) 



APPENDIX B: CALCULATION OF THE 
NORMAL DIELECTRIC FUNCTION AT 
EQUILIBRIUM 



These results are in agreement with |3l| but with the dif- 
ference that we have included the exchange term which 
doubles the interaction strength. 



The expression H28|l can be put in the form of an inte- 
gral over t: 



pi(a)+i04 



/C„(q, w) - 1 



inr-KlJ (Bl) 



Using the development: 



feq 

^k 



(B2) 



we can carry out the integration over the momentum k 
to get: 

/C„(q,c^) = 1 + ^ /°dte'("+'0+)*sin(;Hit) 
m Jq 2m 



3/2 



(B3) 



As far as the imaginary part is concerned the integral 
over the time can be done analytically. The result is: 



ImK;„(q,tj) = -3— rln 

p q 



(B4) 



-/3[m- 



Vl - e 2,,^ 



For the real part the expression, (jB3|l needs to be de- 
veloped in series. For low temperature, Re/C„(q, w) ~ 1. 



APPENDIX C: CALCULATION OF THE 
COLLISION TERM 

As said in section 3, the calculation in the SRPA is 
easy since it presents strong resemblance with the plasma 
gas. The extension to GRPA is more complicated be- 
cause now the condensed and normal modes cannot be 
treated on the same foot anymore, as long as the inter- 
action energy per particle is not the same. Nevertheless, 
the method developed by Ichimaru for a classical plasma 
can be adapted straightforwardly. Here we shall give the 
intermediate steps of the calculations. 

Before doing so, few symmetry properties are interest- 
ing: 

(CI) 
(C2) 



5q(k,k') = .9;(k-q,k' + q) 

^k,k,(q'^))* = ^k,k,(-q,-^) 



The first property can be checked from the definition 
(|36|l and the third one is straightforward. Using the first 
property and (|S7|l . the coUision term defined in Ea. l|l^ 
can be expressed in terms of the imaginary part of the 
correlation function: 

■I ^ 47ra, 



Ck [rik' ; kg 



Let us define: 



^ mV 

q,k' 



Im(5q(k,k')-ffq(k-q,k'))(C3) 



A'(q,c.)=^(l 



(nk + l)nk+q 



k.q 



2m 



(C4) 
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Aq{ci uj) = Hit ~^ ^)^k+q {C5) carrying out the sum over ki and k'^ and taking the ther- 

' UJ + iO+ — ^ — ^ modynamic limit, after a straightforward but lengthly 



With these definitions together with ljA5|l . ljA6|l and 



calculation we find the following expressions for H42|l: 



■ 



duj 



q,k 

Stto , , 



2ni t 



k.q 



31 

2m 



(rik + l)n 



k+q 



/C*(q,c^) 



- 1 



■{l-JC*{q,u))A'{q,u;)-A'*{q,u;) ^ Ao{q,a;)(l - /C*(q,c^)) - A*(q,c^) 



|/C(q,c.) 



|/C(q,c.) 



(C6) 



Sttc 
ml/ 



duj 

q,k' q 

"(l-/C*(q,c.))A'(q,c.)-A'*(q,c.) 



I ?7i 2m 



(nk, + l)nk3+q 



/C*(q,tj) 



- 1 



IC{q,uj)IC*{q,uj) 



_ /C„(q, c^) Ao(q, c^)(l - £*(q, c^)) - ^^(q, m) 
2 |/C(q,c.)|2 



(C7) 



Equivalently X]qk'5q(^s — qjk') is obtained by substi- 
tuting kg by kg — q in (|C7p . The domain of integration 
over w can be extended in the complex plane. Since the 
various reverse dielectric functions in IjASjl converge to 
unity when |cj| oo, the integrand goes to zero in this 
limit faster than l/|w|. As a consequence, the integral 
over to can be carried out by closing the contour either in 
the upper half plane or in the lower half plane. Since, we 
assume that the Bose gas is in a stable regime, the poles 
of the dielectric functions have their imaginary part in 



the lower half of the complex plane, while the complex 
conjugate of these functions have their imaginary part in 
the upper half plane. Thus in ljC6|l and in (jC7|) . contri- 
butions in which poles only lie in one of these half planes 
are canceled since, in that case, the contour of integration 
can be chosen in such a way that no pole is surrounded. 
On the other hand, in order to simplify the Eqs. I|C6|) and 
(|C7|) . we can also add arbitrary contributions in which 
poles lie only in one of these half planes. Taking into 
account these considerations, we arrive at: 



q,K q 



duj 



1 



k.q 



2m 



(rik + l)7ik+q2zIm/C(q, uo) 



Stto 

TTl^lk - rik+qj 

mV 



2ihaA' (q, uj) 2iIm^o (q, t^) 
|/C(q,c.)|2 + |^(q,c.)|2 



(C8) 



E5q(ks,k') 

q,k' 



E 

q 

87ra 



duj 



1 



("ks + l)nks+q2«Im 



(rik, - rtk,+q) 



2m 

2ilmA'{q,u) 
ICiq,Lo)JC*{q,iu) 



IC{q,Lo) 



(1 - ^^^)2zlm 



/C*(q,tj)/C(q,cj 
AQ{q,uj 



\ICiq,coW 



(C9) 



These expressions can be further simplified if we neglect 
some irrelevant infinitesimal terms proportional to 0+. 
Up to these infinitesimal term and using the change of 
variable k' = k -I- q and the formula Im[l/(x -|- iO+)] = 
7r(5(x), we notice that for the two different limits n-i^^/V 



or ni^^/V finite: 

lmlC{q,uj) — ImA;;„(q, cj) 

^ —^22^<'~ci-<')^S{^ - + -^) CIO 

iTiv ^ fn 2to 
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/ IC{q,uj) \ ^ Im/C(q,u;) 



V|/C(q,c.)|2; |/C(q,c.)/C*(q,c.)| 

Finally, it remains to take the imaginary part of these 
expressions, integrate them over lo and plug the results 
into (|C3p . The integration is easy since it involves only 

I 



delta functions. After rearranging terms we notice that, 
in Ea. HC9|l . the term proportional to 1 — A;;„(q, a;)/2 will 
not contribute. Note that, concerning the term (7q(k — 
q, k'), we must carry out the change of variable q — q 
and k' — q — + k'. In this way, we obtain Ea. lB^ and 
Eq.(0. 

Using (jA5|l . a similar reasoning allows also to find the 
expression: 



K,K 



duj 
27ri' 



:2«Im 



- 1 Aoiq,uj) 



(C13) 
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